Emerging Bosons with Three-Body Interactions from Spin-1 Atoms in Optical Lattices 
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We study two many-body systems of bosons interacting via an infinite three-body contact repul- 
sion in a lattice: a pairs quasi-condensate induced by correlated hopping and the discrete version 
of the Pfaffian wavefunction. We propose to experimentally realise systems characterized by such 
interaction by means of a proper spin-1 lattice Hamiltonian: spin degrees of freedom are locally 
mapped into occupation numbers of emerging bosons, in a fashion similar to spin-1/2 and hardcore 
bosons. Such a system can be realized with ultracold spin-1 atoms in a Mott Insulator with filling- 
factor one. The high versatility of these setups allows us to engineer spin-hopping operators breaking 
the SU(2) symmetry, as needed to approximate interesting bosonic Hamiltonians with three-body 
hardcore constraint. For this purpose we combine bichromatic spin-independent superlattices and 
Raman transitions to induce a different hopping rate for each spin orientation. Finally, we illustrate 
how our setup could be used to experimentally realize the first setup, i.e. the transition to a pairs 
quasi-condensed phase of the emerging bosons. We also report on a route towards the realization of 
a discrete bosonic Pfaffian wavefunction and list some open problems to reach this goal. 
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I. INTRODUCTION 

Ultracold atoms and trapped ions offer unprecedented 
possibilities to realize, control and observe quantum 
many-body phenomena [TH3]. For these reasons, in the 
last years they have been frequently employed as quan- 
tum simulators, i.e. controllable laboratory setups mim- 
icking other interesting but not easily accessible systems 
described by the same mathematical model. Simulation 
targets come from diverse research fields, as condensed- 
matter or even high-energy physics. A paradigmatic 
quantum simulator is a system of bosonic atoms in an 
optical lattice, which provides a practically ideal realiza- 
tion of the Bose-Hubbard model @HZ]. On one hand, 
many of these studies aim at the simulation of systems 
which are otherwise difficult to treat numerically, such as 
the Fermi-Hubbard model [5J H] or many-body frustrated 
models 3 J . On the other hand, a parallel goal of quantum 
simulation is the experimental study of "blackboard" the- 
oretical models, such as the Ising model [ID], the Tonks- 
Girardeau gas |11| . or the one-dimensional Dirac equa- 
tion P2]. 

Within the last context, the possibility of simulating 
many-body systems characterized by interactions involv- 
ing more than two particles is of great interest. For 
example, three- and four-body interactions are known 
to be the essential ingredients of lattice gauge theo- 
ries [13] ■ Such theories and their related lattice models 
play important roles in the context of novel exotic quan- 
tum phases, of the breakdown of the Landau-Ginzburg 
scenario and of the confinement-deconfinement transi- 
tion [T3HIS]- Besides that, local many-body interactions 
are essential also in paradigmatic spin-models exhibit- 



*Electronic address: leonardo.mazza@mpq.mpg.de 



ing topological order [T7l419| . Moreover, in the presence 
of external magnetic fields, they lead to various exotic 
fractional quantum Hall states [2UJ [5T] ■ A celebrated ex- 
ample is the Pfaffian wavefunction [55] [53J , which arises 
in bosonic systems at filling factor v — 1 and exhibits 
topological order. Its quasi-excitations are non-Abelian 
anyons, i.e. the exchange of these quasiparticles is as- 
sociated with non-commuting trasformation of the sys- 
tem [24]. Furthermore, Pfaffian states can be found 
within other frameworks, such as p-wave superconduc- 
tivity, where the excitations correspond to zero-energy 
Majorana fermions [55J [55] [55] , or one-dimensional sys- 
tems [27]. This motivates the need for experimentally 
feasible proposals realising three- or many-body interac- 
tions. 

Recently, cold atom theorists have been developing 
several approaches to achieve this goal. The early propos- 
als employed higher-order super-exchange interactions on 
triangular and kagome lattices [55]. Alas, the tempera- 
tures required are even lower than those necessary for 
quantum magnetism in Mott phases, a demanding task 
to which a lot of experimentalists are still committed. 
Super-exchange interactions of the second order involv- 
ing Raman transitions between atoms and molecules in 
square lattices have been proposed to realize an effec- 
tive ring-exchange Hamiltonian for bosons [5J5]. A com- 
pletely different approach has been proposed by Biichlcr 
et al. [30] who suggested to use polar molecules with a 
setup inhibiting two-body interactions. Very recently it 
has been suggested to use the dissipative dynamics in 
presence of a large three-body loss rate in order to en- 
force on the system an effective three-body hardcore con- 
straint [3T] [35] . The idea is very reminiscent of a scheme 
experimentally realized by Syassen et al. |33j to induce 
strong correlations in molecular gases. The dissipative 
scheme has been combined with the rotation of the trap 
to induce artificial magnetic fields and achieve the for- 
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FIG. 1: Sketch of the proposed mapping. A Mott Insulator 
with filling 1 whose atoms have three relevant degrees of free- 
dom can simulate a system of bosons with an infinite contact 
repulsion via the mapping W. The three relevant degrees of 
freedom can be identified with a F = 1 hyperfine manifold. 



mation of the Pfaffian wavefunction [32] . Finally, it has 
been theoretically shown that a perturbative treatment 
of vibrational bands in optical lattices induces effective 
many-body terms in the Hamiltonian, which have also 
been experimentally observed in the time-evolution of the 
quantum phase of the system [331 ES] ■ 

The purpose of this paper is twofold. The first one 
is the investigation of two systems of bosons interacting 
via infinite three-body repulsion in a lattice. We show 
that in one dimension the presence of correlated hopping 
in such models leads to an exotic quantum phase, char- 
acterised by quasi-condensation of pairs of particles. In 
the presence of magnetic fields in 2D, the ground state 
of such model corresponds to the lattice version of the 
bosonic Pfaffian wavefunction. In particular, we discuss 
the stability of the topological properties of such wave- 
function in a non-dilute limit with the magnetic length 
comparable to the lattice spacing. 

The second purpose is to illustrate how spin-1 atoms in 
a Mott Insulator (MI) with filling one offer a possibility to 
tailor three-body interactions in a setup combinable with 
artificial gauge fields. The idea is based on mapping the 
internal states of real spin-1 bosons on the lattice into 
occupation numbers of emerging bosons, similarly to the 
correspondence between spin-1/2 particles and emerging 
hardcore bosons. An easy generalization to higher-spin 
atoms can open the route to the simulation of four-, five- 
and many-body contact infinite repulsions. The first pro- 
posed model, which exhibits the quasi-condensation of 
pairs, can be carried on rather simply in such a frame- 
work. On the other side, the realisation of the lattice 
Pfaffian wavefunction seems to be slightly outside the 
class of models accessible with our proposal. 

The experimental system discussed in this article 
greatly exploits the use of bichromatic superlattices [36] 
in order to individually tailor the hopping rates of the 
different spin species. This scheme is rather general and 



could allow the realization of very general hopping op- 
erators, including even spin-flipping terms; therefore, it 
could find applications beyond the problem of many- 
body interactions, such as in fermionic systems or in non- 
Abelian gauge theories [57] , 

The article is organised as follows: we start in Sec. [Il] 
describing how to realise three-body interacting bosons 
with atomic setups characterized by three relevant de- 
grees of freedom. Subsequently, in Sec. 
sider the explicit case of atoms with F 
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we con- 
1 hyper- 
fine ground manifold and we derive the related super- 
exchange Hamiltonian. We then move to the problem 
of having a complete external access to all the relevant 
parameters of the real spin-1 system and in Sec. |IV| we 
describe how superlattices and Raman transitions can be 
used to achieve this goal. The next two sections are de- 
voted to the discussion of systems characterized by three- 
body interactions and they possible implementation with 
our proposal. In Sec. W\ we present a phase character- 
ized by quasi-long-range order of pairs. 
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we 



In Sec 

discuss the problems faced in trying to engineer the Pfaf- 
fian wavefunction. Finally, in Sec. |VH| our conclusions 
are presented. 



II. THE MAPPING 

As a starting point for our work, we recall here that 
particles interacting via a three-body infinite repulsion 
effectively undergo a three-body hardcore constraint, i.e. 
there can not be more than two particles at a time in 
the same place. Therefore, in the presence of a spatially 
discrete setup, the description of local (on-site) degrees of 
freedom is captured by a finite Hilbert space of dimension 
three: 
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(1) 



Since three-body elastic interactions are rather weak 
in nature, we have to simulate such an ideal system by 
means of a proper realistic discrete setup characterized by 
a local Hilbert space W 1 "^ of dimension three. A unitary 



mapping between the local Hilbert spaces (FigjlJ 
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allows us to relate the real experimental dynamics to the 
dynamics of the emerging bosons characterized by an in- 
finite three-body interaction. If the first is described by 
some effective Hamiltonian H^ al , the correspondent i^hb 
would then be 
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3hb 



W® L ff r f al wt® L 



(3) 



and tuning experimental parameters in -H^al P erm its in 
principle the investigation of a variety of "blackboard" 

-ff3hb- 

In the context of ultracold atoms in optical lattices, 
which offer an unprecedented control on discrete struc- 
tures, it is then quite natural to consider a Mott Insulator 



(MI) with filling 1 and three internal degrees of freedom. 
In this article we focus our attention on an hyperfine 
manifold F — 1, as exhibited by 87 Rb or 23 Na, that have 
already been successfully loaded into an optical lattice 
and cooled into a MI state without freezing the spin- 
dynamics [38 . If the experimental apparatus is tuned 
to a not too deep MI, small quantum fluctuations in the 
atomic position give rise to a non-trivial dynamics, usu- 
ally addressed as superexchange effects. This dynamics 
can be characterized by an effective Hamiltonian -HJeai 
obtained through a second-order perturbative expansion 
of the kinetic term of the Bose-Hubbard Hamiltonian de- 
scribing the real atoms on the lattice [33]. The next 
Section is indeed devoted to a detailed derivation of such 
effective theory. However, we underline that this is not 
the only possible choice and other internal degrees of free- 
dom could have been chosen; as proposed in Ref. [40J for 
a different purpose, these three local degrees of freedom 
could even correspond to different scalar atomic species. 



with m ia 
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b\ a F a fibifi is the total spin on the site i. The first term 
represents the kinetic energy (hopping) and the A Q 's 
represent the energy offset of each of the three states 
(A = 0); the last two terms describe the two-body in- 
teraction on a same site that, due to the spin nature of 
atoms, is characterized by two s-wave scattering lengths, 
a for the S tot = channel and a 2 for the S tot = 2 one, 
with the ratio U 2 /U = (a 2 — a )/(a 2 + 2a ) [31] . 

Hamiltonian Q preserves the total magnetization 
M = J2i &i — lJi a i °f tne sample, allowing us to work 
in a convenient block-diagonal representation. Moreover, 
each energy offset A Q plays the role of a chemical po- 
tential for the atomic species a. As a consequence, in 
absence of spin-flipping interactions (n a conserved), the 
A a would play no role in the dynamics at fixed magne- 
tization. In our case the situation is complicated by the 
presence, in the atomic Hamiltonian, of terms which flip 
the atomic spin: 



III. SPIN 1 ATOMS 

The focus of this Section consists in deriving the effec- 
tive Hamiltonian describing magnetic degrees of freedom 
inside the Mott Insulator and in investigating whether 
interesting 7?3hb Hamiltonians can be effectively mim- 
icked via the proposed mapping. The proposed anal- 
ysis is generally valid for systems of dimension one, 
two and three. In order to describe the experimen- 
tal system of ultracold spin-1 atoms in an optical lat- 
tice, we rely on the standard Bose-Hubbard Hamilto- 
nian |41j (a = { — ,0,+} runs over the three spin states 
{\m F = -1), \m P = 0), \m F = +1)}) : 

#roai = ^ [-t a bl a bj, a + H.c] + ^ A Q n l:Q + 

<z,_7>,ck i,ct 

+ ^5>(n,-l) + ^X;(#-2»,) (4) 
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The A Q would still not play any role in the dynamics if 
the following relation holds: 2A = A + + A_. In pres- 
ence of an external magnetic field, only the linear Zeeman 
shift satisfies this requirement, whereas the quadratic one 
does not. In this case, the relevant dynamical quantity 



is S = A , 



2 A , which quantifies deviations from 



the linear splitting regime. It is experimentally possi- 
ble to control small values of S dressing the system with 
microwave fields [42] . 

In case interaction energies are larger than the hopping 
rates (Uq + U 2 , Uo — 2U 2 >■ \t a \), the system is in a MI 
phase and with a second-order perturbative expansion of 
the kinetic term we compute the super-exchange Hamil- 
tonian i^rcaii wno se link-expression reads as follows: 
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A = 



U 2 (U +U 2 ) 



U2JUa + U2) 



(Uo + U 2 )(U -2U 2 )+SU ^ (U„ + U 2 )(U -2U 2 )+S(U 2 -U ) 
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(U + U 2 )(U„-2U 2 )+5(U 2 -U Q ) 
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Let us now consider a simple class of mappings Ws^x 
which will be used in the following, characterized only by 
a simple phase freedom: 
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These mappings are characterized by the property that 
the magnetization of the spin insulator is directly 
mapped into the density of the three-hardcore bosons. 
Therefore, since the Hamiltonian in Eq. (|6| contains only 
off-diagonal terms which preserve the total magnetiza- 
tion, we automatically gain the possibility of studying 
hardcore bosons setups at fixed density. 

As far as the interaction strengths are concerned, we 
report that the scattering lengths ao and a^ have very 
similar values both in 87 Rb and 23 Na. This means that 
the spin-dependent part of the interaction is in natural 
setups almost negligible, as it is clearly stated once the 
ratio U 2 /U a is calculated, respectively —0.005 and 0.04 
for the two atoms [351 H31 SI] ■ This could in principle be 
a weak point of our proposal, because four non-diagonal 
matrix element in H^eal are proportional to L^. In the 
next sections we will show how to cope with this problem, 
and namely that even in the case in which it is not pos- 
sible to externally tune the ratio U2/U0 to higher values 
it is possible to experimentally observe some interesting 
physics. 

Before concluding this section, we would like to stress 
that the off-diagonal matrix elements present in Eq. ([6]) 
correspond, via the mappings W{ v y in Eq. (Trl), to 
terms describing the hopping of three-hardcore bosons. 
In particular, the matrix elements of the second sub- 
/superdiagonal correspond in i?3hb to one-particle hop- 
ping terms 

|0)|l)o|l)|0), |l)|l)o|2)|0), |1)|2)**|2)|1); (8) 

whereas those of the fourth sub-/superdiagonal corre- 
spond to a two-particle hopping 



|0)|2)*+|2)|0). 



(9) 



The relative importance of these two contributions can 
be modified just by tuning |t |, a factor which multiplies 
the second sub-/superdiagonal and is not present in the 
fourth one. 

In the next sections we will try to use of the Hamilto- 
nian in Eq. |6| to study interesting phenomena related 
to the presence of a three-hardcore constraint: we will 
focus on the realization of a quasi-condensate of pairs of 
particles and discuss the problems in realizing the two- 
dimensional Pfaffian state. Before doing that, however, 
we will describe a setup providing us an high external 
control on the hopping parameters of the system, an in- 
escapable requirement if we want to be able to realize the 
desired Hamiltonians. 




FIG. 2 
V = 20E, 



Plot of the superlattice potential in Eq. ( 10 1 with 
where E T — h 2 k 2 /(2m). We show 



and A = 1.0, 

the exact wannier wavefunctions of the first and third band 
of the lattice, corresponding to the two bound states we will 
use in our proposal. 



IV. THE EXPERIMENTAL SETUP 

We now describe an optical lattice setup, consisting 
of a bichromatic spin-independent potential dressed with 
suitable optical transitions, that once loaded with spin-1 
atoms allow us to experimentally tune the spin-1 Hamil- 
tonian in Eq. Q . Standard spin- independent optical lat- 
tices are indeed not suited for our purposes since they 
preserve SU(2) symmetry, i.e. they naturally give rise to 
perfectly equal hopping rates t a for all the spin orien- 
tations. On the other hand, spin-dependent optical lat- 
tices [H3 SS] suffer heating problems and short lifetimes 
and are thus as well not optimal. At a contrast, our 
setup allows for tuning independently the three hopping 
rates t a without affecting the "traditional" quite long 
lifetimes. In this Section we present a rather qualitative 
description, leaving quantitative analyses in Appendix[A"] 
We start addressing the one-dimensional setup; the gen- 
eralization to more dimensions is sketched at the end of 
this Section. 

The alkaline atoms we propose to optically trap, 87 Rb 
or 23 Na, are characterized by a nuclear spin I = 3/2 
and thus by a ground state splitted into two hypcrfinc 
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TABLE I: Numerical values of the parameters of the super- 
lattice in Eq. flu) with Vo = 20E r and A = 1.0. We first list 



the energies of the lowest localized wannier functions of main 
and secondary minima and their energy difference; we then 
calculate the parameters of the main lattice, i.e. the hopping 
rate J ma in-main and the interaction energy Uon 
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FIG. 3: Sketch of the scheme we propose to induce hopping 
between the levels of the F = 1 manifold, i.e. the adiabatic 
elimination of one F = 2 state trapped in the intermediate 
minimum (red non-dashed arrows). Because of orthogonal- 
ity properties of wannier functions, the coupling cannot be 
realised with microwave fields. Optical Raman transitions 
through an excited state carry non-negligible momentum and 
can therefore be a solution. In Appendix |A| we also discuss 
the effect of spurious couplings as those depicted with orange 
dashed arrows. 



manifolds F — 1 and F — 2. If the laser originating 
the lattice is detuned enough from the excited levels, 
the optical potential is insensitive to atomic spin prop- 
erties and the eight ground hyperfine levels experience 
the same dipole potential [47]. By combining two lasers 
with exactly commensurate frequencies 1:2 (obtainable 
by standard frequency-doublers) it is possible to create a 
bichromatic superlattice 

V(x) = -V [cos 2 (fca;) + A cos 2 (2kx)] ; V ,X>0, 

( 10 ) 
where Vo is proportional to the overall light intensity and 

A describes the relative strength of the two lasers. Tuning 
these two parameters, it is possible to create a configura- 
tion with main (deep) and secondary (shallow) minima, 
as depicted in Fig. [21 In particular, it is experimentally 
possible to cool a MI with atoms only in the main min- 
ima just superimposing the second lattice after the cool- 
ing has already been done. We propose to store all the 
physical information into the main minima of the F = 1 
manifold whereas the other one will only provide auxil- 
iary levels. Indeed, the idea is to use the levels trapped in 
the secondary minima of the F = 2 manifold to provide 
intermediate "bus" states to be adiabatically eliminated, 
as pictorially explained in Fig. [3] 

The bound states energies and the corresponding 
localized Wannier wavefunctions were calculated solving 
the single-particle Schrodinger equation with periodic 
boundary condition and manipulating the periodic 
solutions [48] . The overlap between them is exactly zero 
because they belong to different bands of the lattice 
spectrum, making it impossible to implement couplings 
via microwave fields that carry negligible momentum. 
Our suggestion is to engineer them via optical Raman 
transitions, adiabatically eliminating a far excited state 
like the one employed to create the trapping optical 



FIG. 4: Splittings of the levels of the F = 1 and F = 2 
hyperfine manifolds in 87 Rb due to an external magnetic field. 
The splitting between the two manifolds is not in scale. Red 
arrows describe the effective couplings we want to engineer 
via Raman transitions, <5s and ils are the effective parameters 
describing these transitions. 



dipole potential and consequently transfering a momen- 
tum comparable to the lattice inverse spacing. Provided 
a non-null intersection of the effective support for the 
Wannier wavefunctions, this scheme permits to have 
couplings like those of Fig. [3] We note that, because 
of the coherence properties of laser light, it is possible 
to give a complex phase to the effective hopping and 
therefore to combine this setup with current proposals 
of artificial abelian and non-abelian gauge fields [4"5hTjT] . 
The results presented in this rticle were obtained by 
fixing the parameters Vq — 20E r and A = 1.0, where 
E r = h 2 k 2 /(2m) is the recoil energy of the long-wave- 
lattice; some important parameters of this lattice are 
listed in Table [Q 

We suggest to use the levels of the \F — 2, mp — m) as 
ancillas for the \F = l,mp = m) states. The possibility 
of individually addressing each of the three transitions 
|l,m) -H- 1 2, to) is granted by an external magnetic field, 
which splits the levels according to the formulas: 

Ei, m = -Qf hbB m + A\' m (iJ,BB) 

E 2 , m = &hf +g F fJ, B Bm + A {2 ) n (^ B B) (11) 

where \lb is the Bohr magneton, gp is the hyper- 
fine Lande Factor and Ahf the hyperfine splitting (see 
Fig. Uk. A)/ m (/z_B-B) denotes the second order Zeeman 
shift proportional to ((XbB) 2 of the magnetic sublevel 
mp = m of the spin-i* 1 manifold, which in 87 Rb has 
opposite signs for the two F — 1 and F — 2 mani- 
folds. The values for 87 Rb are fi B gF = 0.7MHz/G and 
Ahf — 6.8GHz. The transitions can be therefore de- 
tuned to a regime in which spurious effects can be safely 
neglected. For example, if we consider 87 Rb, even weak 
magnetic fields of 10 — 100 G can detune the three transi- 
tions of 15— 150 MHz. A more quantitative and detailed 
discussion of these ideas is given in Appendix [Aj where 
analytical and numerical arguments are provided in sup- 
port of our proposal. In particular, we discuss the effects 



of spurious couplings between other localised states of the 
two manifolds, represented in Fig. [3] by orange dashed 
lines, which prove not to affect the efficient population 
transfer between neighbouring sites. 

As far as the generalization to more dimensions is con- 
cerned, we simply propose to apply the laser configura- 
tion originating the potential in Eq. ( f 1 also in the other 
directions, labelled by i: 



V(x) 



-V o y~] [cos 2 (kxi) + \ cos 2 (2kxi)] ; V ,\ > 0; 



(12) 
An increased number of dimensions makes the structure 
of the relative minima more complicated, however it is 
still possible to recognize a square or cubic geometry 
of main minima and an "auxiliary" lattice of secondary 
minima trapped in the middle of links between the main 
ones. Other higher energy minima appear at the centers 
of faces and cubes, but they can be neglected due to their 
even higher energy offset. Similarly to what done in one 
dimension, the system must be then dressed by optical 
lasers driving the nearest neighbor hopping. 

Finally we notice that the other parameters appearing 
in Eq. Q can be experimentally tuned with technologies 
standardly used in current optical lattice setups. The 
interaction strengths Uq and C/2 can be indeed modified 
via Feshbach resonances of the scattering lengths a,Q and 
0,2 [21 [52] whereas the energy offsets A a can be varied 
using the Zeeman effect or dressing the levels with far- 
detuned microwave fields [42] . 

Before concluding this Section, we would like to remark 
that we are suggesting to observe global magnetic proper- 
ties arising from the super-exchange effects of a spinorial 
Mott Insulator. Currently a lot of experimental efforts 
are devoted to this task, mainly within the context of the 
Fermi-Hubbard model, the challenge being represented 
by severe temperature and entropy requirements [53H55J . 
Nonetheless, the huge quest going on makes us believe 
that such phenomena will be experimentally achievable 
in the next-future. In the same spirit we intend also 
the superlattice setup, which at a first glance could seem 
rather intricatcd. 



and of a correlated two-particle term Q in the emerg- 
ing Hamiltonian -ff 3 hb of Eq. (|6| . As already pointed out 
at the end of Sec. |HI[ the relative strength of these two 
terms can be varied by just tuning the strength of \t \: 
with the help of the setup discussed in Sec. IV it becomes 
feasible. We explicitly discuss the realistic case of Eq. ffib 
where many spurious terms emerge but do not prevent 
PQC to be observed. 

At zero temperature (T = 0), very general theorems 
state that no long-range order can arise in ID and Bose- 
Einstcin condensation is consequently ruled out; any- 
way, the presence of algebraical decays in the density 
matrix allows one to introduce the concept of quasi- 
long-range order and quasi-condensation [56] . Beside the 
usual "atomic" quasi-condensate (AQC), characterized 
by quasi-long-range order of the one-particle density ma- 
trix (aja,j + A), it is possible to speak of "pairs" quasi- 



t 2 „ 2 



i+A/ 



condensate when the two-body density matrix (a] a 
still exhibits quasi-long-range order despite the exponen- 
tial suppression of one-particle correlations. 

Even if mean-field calculations support the conjec- 
ture that PQC induced by correlated hopping is not 
a low-dimensionality phenomenon, the simplest setup 
for both experimental and numerical purposes is of- 
fered by a one-dimensional lattice. On one side, a re- 
duced number of Raman pairs of beams is required; on 
the other, numerical simulations can be routinely done 
with the help of Density Matrix Renormalization Group 
(DMRG) [57] [55] ■ In this paper we made use of an open- 
source code (www.dmrg.it) with open boundary condi- 
tions (OBC) but in principle also parabolic external po- 
tentials - closer to current experimental setups - lie within 
the possibilities of the method [58- . 

Before starting our discussion, we report that the same 
paired phase has been also at the focus of Ref. [3T] where 
the transition was driven via two-body attractive inter- 
actions and the three-body hardcore constraint was ef- 
fectively induced by strong dissipation channel. 



A. PQC Induced by Correlated Hopping 



V. PAIRS QUASI-CONDENSATION 

In this Section we study the first interesting model 
characterized by three-body infinite repulsion: we show 
that a dominating correlated hopping can drive a tran- 
sition to a quasi-condensate of pairs (PQC), without the 
need of any two-body attraction. Moreover, we show 
that substituting the three-body interaction with a two- 
body one the system becomes unstable towards collapse: 
this strictly links the PQC to the stabilizing effect of the 
three-body repulsion. 

The experimental realisation of such a phase with the 
help of our setup is possible. Besides the 3-body hardcore 
constraint, we already highlighted in Sec. |III| the simul- 
taneous presence of the usual single-particle hopping ([8]) 



In order to understand whether a system with domi- 
nating correlated hopping undergoes a phase transition 
to PQC, we start studying this simple "blackboard" 
Hamiltonian on a one-dimensional setup with L sites: 



H = - JJ2 ajai+i -KJ2 a ^ a2 



t\3 



i+l 



+ H.c- (a]Y = 0. 

(13) 
The correspondent phase diagram, plotted in Fig. [5] dis- 
plays a large region characterized by an exponential de- 
cay of the particle-particle correlator (ajaj) and by an 
algebraic decay of the the pair-pair one (a] a 2 ). A mean- 
field analysis via Gutzwiller ansatz shows the presence of 
the same phase transition and thus supports the robust- 
ness of the effect even in larger dimensions. It is also 
possible to Fourier transform both density matrices and 
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FIG. 5: Phase diagram of the one-dimensional model in 
Eq. ( |13| ) calculated with a DMRG algorithm. Two phases 
appear, characterized respectively by quasi-long-range order 
of the one-body density matrix (AQC) and by exponential 
decay of the one-body-density matrix and quasi-long-range 
order of the two body one (PQC). The red dashed line is the 
mean-field result obtained via Gutzwiller ansatz. 
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FIG. 6: Phase diagram of the system described by Eq. (141 
in which the three-body interactions have been substituted 
by two-body ones. Instead of a PQC a phase appears which 
is unstable towards collapse. At incommensurate filling and 
J = the ground state manifold is largely degenerate and is 
spanned by Fock states with less than double local occupan- 
cies (Mott Glass). 



analyze the finite-size scaling of the population in the 
largest occupied state, i.e. the lowest momentum one. 
An algebraic growth ~ L a is another signature of quasi- 
long-range order, and we find agreement between the two 
benchmarks. In the following, we only look at the decay 
of the correlators. 

We stress that the only presence of correlated hopping 
is not enough to create a PQC and that the stabilizing 
action of a three-body hardcore constraint is for this sake 
crucial. For example, we can study a model in which we 
substitute the constraint with two-body repulsions: 



H = -J 
-K 



J2 K a ' 

i 

J2 [ a ? a2 *+i 



, , i -r H.c 
H.c 



U^mim-i) (14) 



The relative phase diagram is plotted in Fig. [6] for two 
different densities, n = 1.0 and n = 0.75: when K domi- 
nates a phase appears which in the thermodynamic limit 
is unstable towards collapse. The instability is induced 
by the correlated hopping term and has been already ob- 
served in other numerical works [55] ; in Eq. (13) it was 
counterbalanced by the hardcore constraint. The sta- 
bilizing action of strong two-body U drives the system 
outside this instability region, but alas prevents also dou- 
ble occupancies and then the desired PQC. In the limit 
of large J/K the usual Bose-Hubbard physics made of 
AQC and MI only is recovered (not shown in plots). 



B. PQC in Spin-1 Mott Insulators 

Coming back to the experimental setup proposed in 
Sec. ~ 



quasi-condensate. In order to do that, we use a proper 
mapping (TtJ) Wpqc with all the phases set to a same 
value ip a — 0. The diagonal terms in Eq. (|6| describe 
additional two-body and nearest-neighbors interactions 
that we cannot get rid of; therefore the Hamiltonian in 
Eq. ( 13 1 cannot be easily recast and additional numer- 



ical simulations are needed to characterize our approxi- 
mation. 

All the next plots share the same value of t±/Uo = 0.1 
and of S = 0, i.e. no relevant quadratic corrections to 
the linear Zeeman splitting. With respect to the ratio 
U 2 /U , we studied the cases of 87 Rb (-0.005) and 23 Na 
(+0.04) as well as the value U 2 /U Q = -0.04. The phase 
diagrams, shown in Fig.[7J have been studied varying the 
total density of three-hardcore bosons, i.e. the magneti- 
zation of MI, and the ratio between single-particle and 
correlated hopping: 



J 
K 



2t t 



+1 



U0 + U2 



U + U 2 y/2t+it- 



10V2 



u 



(15) 



For the values of the two alkalis the system clearly 
exhibits a PQC phase, even if no clear signature of a AQC 
phase has not been found. Instead an inhomogeneous 
phase appears characterized by phase separation between 
fillings — 1 or 1 — 2. However, a slight tuning of U2/U0 
to —0.04 could already help the system recover all the 
interesting physics of the Hamiltonian in Eq. (13), 



III we investigate whether it supports the pairs 



shown by the phase diagram in Fig. [7J calculated for such 
value. In this case the matrix elements of the second sub- 
/superdiagonal arc still far from being proportional to 
the exact values {1, y/2, V% 2}. The phase transition is 
observable because the finely tuned values of the matrix 
elements are not important, as far as they roughly share 
the same order of magnitude and this last value is varied 
with respect to the correlated hopping parameter. 
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FIG. 7: Phase diagrams of the Hamiltonian Hshbc real- 
ized with our proposal which approximates the model in 
Eq. ( |13[ ). The two plots are drawn for the realistic values 
of U2JU0 = -0.005 ( 87 Rb, left) and U 2 /U = -0.04 (right); 
the case U2/U0 = 0.04 ( 23 Na) is not shown since it is quali- 
tatively equivalent to the case of 87 Rb. On the left, even if a 
PQC phase appears, there are no signatures of AQC. Instead, 
at density n = 1.0 we find a MI whereas at n 7^ 1.0 an inho- 
mogeneous phase appears (phase separation). On the right, 
the phase diagram shows that even a small tuning of U2/U0 
from the atomic values let the AQC phase arise. 



The density profile and the spatial decays of the 
particle-particle and pair-pair correlators for the system 
at n = 1.125 and J/K — 1.54 (AQC) are shown in 
Fig. [8j whereas in Fig. [9] the plots refers to n = 1.125 
and J/K = 1.43 (PQC). The plots show quite clearly the 
presence of a region in the phase space where the pair- 
pair correlator exhibits quasi-long-range order whereas 
the particle-particle one is exponentially suppressed. The 
accurate definition of the phase border, requiring numer- 
ics on larger systems and finite size scalings, lies beyond 
the purposes of this article, which only aims in determin- 
ing the presence of a PQC phase. However, simulations 
for large systems up to 240 sites show that the PQC phase 
is indeed stable and is not a finite-size effect. 

For a more accurate theoretical treatment of the prop- 
erties of the transition from QC to PQC we refer the 
interested reader to Ref. [BU]. The applicability of such 
theoretical methods to our system must not be taken 
for granted because in our setup the phase transition is 
induced by correlated hopping whereas in the referred 
papers by two-body attractive interactions. 



C. Experimental Observation 

Since one of the key features of a quantum simula- 
tor is the possibility of observing the quantum state that 
have been realised, we now discuss how the two AQC and 
PQC phases could be detected with our setup. For this 
purpose, we translate into the spin language the corre- 
lators (a\aj) and (0] a 2 ) which we used to identify the 
two phases. Once the mapping Wpqc is considered, the 
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FIG. 8: Plot of the density profile (n(x)) of the ground state 
of the system for U 2 = -0.04 U , n = 1.125 and J/K = 1.54. 
The two insets show the exponential decay of the particle- 
particle and pair-pair correlators (log-log plots). These data 
allow us to identify the phase as a AQC. 



three-body hardcore operators can be rewritten as fol- 
lows: 



y/2- 1 

^2~ 



*-\: £>« 



(16) 



The pair-pair correlator assumes therefore a very sim- 
ple expression, whereas the particle-particle one can be 
written as: 



(4 a j) = 



V2-1 

3-2V2 



[(SfSfSr) + (S+S-S*)} 



(S i b t Sj bj ) 



(17) 



Numerical simulations show that the analysis of the de- 
cay of the spin correlators (S^~S~) and (S^~ 2 S~ 2 ) leads 
to definition of the same phase boundary as before. The 
possibility of using spin-spin correlators to identify the 
phases is experimentally relevant, because this is the 
most natural language to analyze the properties of a 
Spin-1 Mott Insulator. Moreover, we think that the ac- 
tual significative experimental efforts in order to develop 
techniques able to resolve the single sites of optical lat- 
tices [5TJ [S5] will make the direct observation of the pro- 
posed correlators possible. 



VI. PFAFFIAN PHYSICS 

In this Section we move to the analysis of a second in- 
teresting many-body system characterized by three-body 
interactions: the Pfaffian wavefunction [23]. This state 
has been proposed in the context of the Quantum Hall Ef- 
fect (QHE) [2U] in order to describe the many-body elec- 
tron liquid at fractional magnetic filling v — 5/2. The in- 
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FIG. 9: Plot of the density profile (n(x)) of the ground state 
of the system for U 2 = -0.04 [To, n = 1.125 and J/K = 1.43. 
The two insets show the exponential decay of the particle- 
particle correlator (log plot) and the algebraic decay of the 
pair-pair correlator (log-log plot). The presence of this last 
quasi-long-range order allows us to identify the phase as a 
PQC. 



terest in this wavefunction resides in the predicted prop- 
erty of supporting non-Abclian quasi-excitations [221 163) . 
Here we deal with the bosonic version of the Pfaf- 
fian state (y — 1) and show that this wavefunction 
can be studied also in a lattice. Combining exact- 
diagonalization numerical approaches and some well- 
known benchmarks to test topological properties, we see 
that even at significative magnetic fields the ground state 
of the system features non-trivial topological hallmarks. 
We then employ these tools to discuss the possibility of 
using a spin-1 MI to realize such a wavefunction, and 
underline some still present drawbacks in the recipe. 



A. Quantum Hall Effect on a Lattice 

We consider a two-dimensional setup with 7Y bosons 
with charge q interacting via purely three-body repulsion 
(no two-body term) in presence of an external uniform 
magnetic field with vector field A. The setup is pierced 
by a number of magnetic fluxes iV$ equal to the number 
of particles N (filling factor v — 1); a typical length 
I = ^Jhc/qB is induced in the system by the magnetic 
field itself. The system is ruled by the following many- 
body Hamiltonian, in which we write the position of the 
particles with complex coordinates z = (x + iy)/£' 



H 



pf 



E 



l A( Zl ) 



2m 



■<*E 

i<j<k 



5(zi - z j )8{z i - z k ). 

(18) 

C3, greater than zero, is the strength of the repulsion. The 
single particle levels are arranged into a collection of de- 
generate manifolds, the Landau Levels (LL), separated 
by a gap twice the cyclotron frequency 2h(qB /mc); as 
long as the chemical potential is smaller than this separa- 



tion, the particles will live only in the lowest LL and will 
be characterized by wavefunctions analytical in z (the 
exponent being the angular momentum). Within this 
framework, the double-delta potential is properly regu- 
larized and the ground state of the Hamiltonian is the 
Pfaffian wavefunction 231: 



*(zi,...z w )ocPf 



1 



n^ ~ *j) 



■Eji^r/2 



i<3 



(19) 

The Slater determinant IL^^ — z j) w °uld prevent the 
coincidence of two or more particles in the same spa- 
tial position; the prefactor Pf(l/(zj — Zj)), the Pfaffian 
(square root of the determinant) of the antisymmetric 
matrix with elements A^j = l/(zi — Zj), enables the su- 
perposition of two bosons but still forbids that of three. 
With this construction, the wavefunction is forced to be 
the lowest angular momentum state in the intersection 
between Lowest Landau Level and kernel of the three- 
body interaction. 

In order to discuss the possibility of simulating the 
Pfaffian state with our proposal, we have first to dis- 
cretize the system. We take into account the presence 
of a three-body interaction with c% — > oo introducing the 
three-hardcore bosons operators a and a^ satisfing a 3 = 
and a^ 3 = 0. The discrete version of a kinetik Hamilto- 
nian with minimal coupling is: 



Hpf_ 



Pf-lat 



-J £ e^alaj+H.c; (a]) 3 = 0. (20) 



<hj> 



As in every discrete U(l) gauge theory, the magnetic field 
coupling to the positional degrees of freedom of the par- 
ticles is represented by a phase ifii.j — 27r/$o f J A ■ dl, 
where $o = hc/q is the quantum of flux. 

We already discussed in Sec. |H] and |IH| how to exper- 
imentally deal with the lattice version of three-hardcore 
bosons; we also showed that our superlattice setup is 
compatible with the general theoretical idea of inducing 
a phase in the hopping with an electro-magnetic running 
wave as in the pioneering proposal by Jaksch et al. [1H] ■ A 
flurry of theoretical proposals and experimental attempts 
to realize an artificial gauge field for neutral atoms have 
been going on for years now [SOI [STJ 151-68 and have 
been further spurred by the breakthrough work by Lin 
et al. [69) that illustrated the experimental realization in 
a BEC. Therefore we think that also this technological 
aspect of our proposal lies within the next-future possi- 
bilities and the only problem we are left with is whether 



it is possible to realize the model Hamiltonian ( 20 1 within 



the framework of Spin-1 Mott Insulators with Raman su- 
perlattice dressings. 



B. Topological Properties as a Benchmark 



Before discussing the simulation of Hamiltonian ( 20 1 , 



we investigate to which extent transposing the physical 
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system onto a discrete lattice modifies the nature of the 
many-body state. The problem arises from the competi- 
tion of two typical lengths, the magnetic one I and the 
lattice constant a. In the small magnetic field limit I 3> a 
(or dilute limit, since the constraint N = iV$ must hold), 
we expect the system to be insensitive to the discrete 
nature of the space. On the other side, an analysis of 
what happens when the magnetic field (and the particle 
density as well) increases is needed to test the robustness 
of a fully discrete version of the Pfaffian wavefunction. 

The characterization of QHE wavefunctions trans- 
posed from continuum systems (usually two-dimensional 
strongly interacting electrons) to discrete optical lattices 
is a problem that has already been faced in the literature 
[70TfT2] . Here we follow the standard approach. We per- 
form an exact diagonalization of the system with periodic 
boundary conditions (PBC). Three marks are used to test 
the genuine Pfaffian nature of the numerical ground state 
(see Table In] for their values in this case) : 

1. the agreement between the degeneracy of the dis- 
crete numerical and continuum analytical ground 
manifolds (the Pfaffian wavefunction has been gen- 
eralised on a torus first in Ref. 23J); 

2. a significative overlap of the discrete numerical 
wavefunctions with the continuum analytical ones; 

3. the agreement between the Chern number 
(CN) [731 El] of the discrete numerical and 
continuum analytical ground manifolds. 

We stress that the three-fold degeneracy of the pfaffian 
ground state is not of topological nature, and is strictly 
connected to the properties of the Jacobi theta functions, 
which are used to generalized on the torus some QHE 
states |!?3"] . 

Chern Numbers probe the topological properties of 
the system testing its sensibility towards the twist of 
the boundary conditions, expressed by two parameters 
{6x,6y) G [0, 2tt) x [0, 2tt). We give here the expression 
of the first CN for the simple case of non-degenerate 
ground state, whereas for more dimensions we refer to 
Refs. [7J3I7J]: 



C = 



^J d9 x d9 v [de^Ay^ez) - d 6y A x (ei,6 2 )] (21) 



where A t = (V(0 X , e y )\^-\^{9 x , 6 y )) and |tf(0 X) v )) is 
the ground state with boundary conditions {9 X , 9 y )- This 
integer quantity is indeed related to the theory of topo- 
logical invariants in the context of the Berry connection. 
CN are increasingly used in condensed-matter theory 
since the discovery that the quantized properties of the 
anomalous QHE resistivity could be studied within such 
framework [75]. We calculate the CN with the method 
provided by Hatsugai |73j which avoids any explicit nu- 
merical differentiation and connects the CN evaluation 
to the number of vortices displayed by a proper auxiliary 
field £l(9 x ,9 y ). The consequent integer character of the 




ft, 



FIG. 10: Plot of the auxiliary field Q,(6 x ,9 y ) for the system 
with Hamiltonian (20 1. The parameters of the system are 



those in Table |H| The three highlighted vortices mean that the 
CN of the system is equal to 3. Vortices must be searched in 
regions where the main chosen gauge has minima: the contour 
lines and the shading highlight such part of the parameter 
space. The definition of the field Q. and the way it can be 
computed are discussed extensively in Ref. [70] [71] [73] , to 
which we refer the interested reader. 



CN constitutes a further reason of its extensive use, since 
it provides a reliable yes-no benchmark more immediate 
than the wavefunction overlap ranging in [0, 1]. 

Within the uncertainty given by working with small 
systems without accessing the thermodynamic limit, we 
can at least affirm that our results are compatible with 
the presence of an incompressible liquid with a degener- 
ate ground state on the torus at £ ~ 0.8 a. Moreover, 
they also present significative signatures that the nature 
of the system should be strictly linked to that of the Pfaf- 
fian state. 



Tentatives Towards the Clean Model 



The previous Section shows that if we were able to im- 



plement Hamiltonian ( 20 1 we would access the intriguing 
physics of Pfaffian with our quantum simulator. The 
Hamiltonian for effective bosons obtained in Eq. ( |6"| h as 



to be then compared with the link version of Eq. ( 20 1 
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where the graphic symbols highlight some terms of 
Eq. ([6| which are not present in ( 20 1 . 
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U 2 = V2/(2V2 + 3)U ~ 0.24J7 <5 = -2\t o \ 2 /{U + U 2 ) 



to = 0.W 



t- = O.ltoe^ t+ = 2il 



TABLE II: Exact diagonalization study on a torus of the 
many-body ground state of the system described by the 
Hamiltonian in Eq |20| The degeneracy and the Chern number 
of the ground manifold in the continuum case are respectively 
3 and 3. As discussed in Appendix [B] the presence of mag- 
netic fields strongly constraints the dimension of the torus to 
be simulated; the next size would be 5 x 5, with an Hilbert 
dimension 110k. 



TABLE III: Set of parameters used together with the mapping 
Wpf to recover the model in Eq. |20| We stress that it is 
possible to give to the phase •& the space dependence which 
characterizes d>i ,. 



Unfortunately we were then not able to find neither a 



One of the problems is related to the presence in <\? 
of the correlated hopping term p)]), which is not com- 
prised by the the QHE model. Therefore, we tried to 
study the model in the regime: |i | » \t+\, |i_|, which 
decreases the relevance of correlated hopping. In this 
case we use a mapping Wpp characterized by the phases: 
{ip- = 0; ip = 0; tp+ = it} and set the various parame- 
ters to the values listed in Table IIIII This sets the second 
sub-/superdiagonal to be approximately proportional to 
{1; \/2; V2; 2}. Unfortunately, this tunes only eight of the 
terms of the diagonal to an approximate same value: the 
central one A is significantly different from the others, 
leading to a completely different model with an effective 
nearest neighbours interaction. Moreover, this method 
has the general disadvantage that the effective hopping 
rate J would be proportional to \t t-\/(Uo + U-z) and 
therefore require temperatures even lower than the pure 
super-exchange effect \t Q \ 2 /{Uq + Uz). 

As an alternative, we abandon the attempt to exactly 
recover the model in Eq. (20) and try instead to realize 



way to get the Hamiltonian in Eq. (20) nor to realize a 



a similar system whose ground state is characterized by 
the same benchmarks of the Pfaffian wavefunction, i.e. 
the same degeneracy on the torus and the same Chern 
number. At low density, the number of global Fock states 
with more than two particles on one link is lower than 
that of the other states. Thus, we expect that matrix 
elements of the link Hamiltonian connecting states with 
more than two particles per link do not play a relevant 
role in the global dynamics; even sensible deviations 
of such terms from the exact values should not change 
too much the properties of the ground state. Hence, we 
investigated sets of parameters which could put all the 
"noise" on such matrix elements. We consider the same 
4x4 system as before at density p = 1/4 and magnetic 
field iV$ = 4, which we can numerically analyze, but 
the next considerations could also be generalized to 
systems with smaller magnetic fields (or more dilute). 
Unfortunately, even this turned out to be impossible. 
We tried to combine a tomographic analysis of the 
Pfaffian wavefunction with the tuning of all the matrix 
elements of the link Hamiltonian connecting states 
with less than three particles. However the numerical 
simulation of these Hamiltonians gave always as result 
non-degenerate ground states characterized by no topo- 
logical properties, i.e. a Chern number equal to zero |74j . 



similar Hamiltonian whose ground state was three-fold 
degenerate and characterized by a Chern Number equal 
to three. Ergo we think that, alas, the Pfaffian wave- 
function cannot be readily implemented with the help 
of a quantum simulator based merely on the ingredients 
described in this work. It might be nonetheless the case 
that, adding to the proposed setup some furher trick or 
ancillary system, it becomes feasible. 



VII. CONCLUSIONS 



In this article we have discussed two many-body ex- 
amples of systems characterized by a three-body infi- 
nite contact repulsion. In the former we studied a one- 
dimensional phase characterized by quasi-long-range or- 
der induced by correlated hopping, whereas in the latter 
we have examinated the stability of a discrete bosonic 
Pfaffian wavefunction in a non-dilute limit. 

Moreover, we suggest to experimentally realise such 
phases with the help of optical lattices and spin-1 atoms. 
Our proposal to simulate three-body infinte repulsion re- 
lies on a local mapping between the dynamics of a spin-1 
MI and that of emerging bosons characterized by such 
interaction. Numerical calculations support the experi- 
mental feasibility of the former setup, whereas the latter 
it seems that further control parameters are still needed. 

A crucial point of this paper is the extensive descrip- 
tion of a bichromatic optical superlattice which could al- 
low the realisation of rather general hopping operators 
for spin gases in optical lattices. We describe the case of 
spin-preserving hopping rates and show that it is possi- 
ble with laser assisted tunneling to break the SU(2) sym- 
metry. Moreover, in a next work we will show that an 
additional staggering of the lattice can in principle give 
access to a more general class of hopping operators oper- 
ators containing even spin- flipping terms. We hope with 
this to open a route towards theoretical studies of inter- 
esting models or towards the experimental realisation of 
exotic spin models. 
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Appendix A: External Control of the Hopping Rate 
via Superlattices 

In this Appendix we give a quantitative analysis of the 
discussion of Sec. |IV| in which we suggested to use super- 
lattices in order to externally and independently control 
the hopping of the three spin species. 

We start discussing the explicit expression of the cou- 
pling realised with an optical Raman transition between 
two different hyperfine levels of the ground state L = 
via elimination of the manifold of excited states L = 1, 
where L is the electronic angular momentum. Atomic 
levels are addressed with the notation \L, a, k), where 
a labels the hyperfine degrees of freedom and k are 
the quantum numbers of the center-of-mass wavefunc- 
tion. The Raman coupling between two states \0ak) 
and \0a k) is: 

fW ;a *(t) = -\ Y. (k'\e- ip2 '*\q)(q\e ipi -*\k) ■ 



|l/3?> 



_ —i(0Jl— U!2)t . 



' c *2a> 13 \\V\\*2 E 2 E 1\\V\\1 ClaP ' e 
1 1 

Elfiq — E 0a k — hwi Eip q — -Efla'A;' — ^1 



(Al) 



where ELak is the energy of the level \L, a, k) and u>i and 
Pi are the energy and momentum of the i-th laser. The 
coupling realised by the z-th laser between the internal 
atomic states \0a) and |l/3) is described by ||/x||i, E t 
and ci a p according to the notation of Ref. [37]. The 
sum over the excited states is limited to the first excited 
manifold because we consider lasers far-detuned from 
higher excited levels. In the case of a spin-independent 
lattice, lasers can be detuned from the first excited 
manifold L = 1 of even some tens of THz: in this case 
the expression in Eq. (Al) can be simplified. Indeed, 



the energy differences at the denominators depend only 
slightly on the internal structure of the levels (they 



can differ at most for some GHz): once E\p q — Eo a 'k> 
is substituted with the 0-th order energy difference 
between excited and ground states AE W , we can write: 



i ({W\t 



*(P2-Pl)-X 



\k) , {k'\e 



'|p- i (P2-Pl)-X 



AEio — hu>i 



+ 



\k) 



A_Ein — fVaJi 



•J2 C 2a>{3 IIHI2 E 2 ElWlAllClap ? 

& 
q q —tut 



— i(LOl—UJ2)t 



(A2) 



In this last expression u = cJi — W2, Sk<k — 
(/c'|e _I< - P2_Pl ' > x |fc) whereas Q a > a comprises all the re- 
maining terms. The very simplified expression for the 
center-of-mass part of the coupling Swk comes from the 
substitution of J2 q \q)(q\ with the identity on the center- 
of-mass Hilbert space. 



Taking advantage of Eq. A2 specified to the setup de- 
scribed in Sec. |IV[ we now discuss the possibility of trans- 
ferring population between two F = 1 neighbouring sites 
via adiabatic elimination of an F = 2 state trapped in 
the middle. 

We consider two states with the same magnetic quan- 
tum number tof, \F — 1, nip) and \F = 2, my), and de- 



velop the "6- levels model" depicted in Fig.[TT] We believe 
this model captures the relevant physics of superlattices 
dressed with one Raman coupling and includes spurious 
couplings between main and main or secondary and sec- 
ondary sites. The coupling between levels trapped at 
different positions, i.e. belonging to different bands of 
the lattice, is possible only because we are transferring 
momentum via the lasers. Since we are working at fixed 
L = and mp, we restrict the previous notation \0ak) 
to the shorter \F, k), the two quantum numbers being the 
hyperfine manifold F — 1,2 and the position where the 
center-of-mass wavefunction is trapped (for the mean- 
ing of k = 1,2,3 see Fig. 11). The model is character- 
ized by only three relevant Sk'ki as depicted in Fig. [XT] 
couplings between neighbouring main sites are negligible. 
The Hamiltonian reads as follows: 

H= d|l,2)(l,2|+A(|2,l)(2,l| + |2,3)(2,3|) + 

+ (A + d)|2,2)(2,2| + 

+ Qe-^l S 1 , 2 (|2,2)(l,l| + |2,2)(l,3|) + 
+ 5 M (|2,1)(1,1| + |2,3)(1,3|) + 
+ S 1 *, 2 (|2,1)(1,2| + |2,3)(1,2|) + 

+ S 2 , 2 |2,2)(1,2|] + H.c. (A3) 

Once we apply the unitary transformation r(t) = 
exp[£d(|l,2)(l,2| + |2,2)(2,2|)t], the three levels |l,fc) 
become degenerate. In case the three inequalities 
\Sijfl\/(5 — d) 4C 1 arc fulfilled, it is possible to use 
second-order perturbation theory in order to develop an 
effective Hamiltonian describing the dynamics within this 
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1.00 



F=2 



F=l 



I* = 1) 



\k = 3) 



FIG. 11: The "6-levels model" used to study the coupling of 
different hyperfine levels with one Raman transition. Levels 
are labelled by two quantum numbers, F and k. Energies are 
not in scale; the orders of magnitude of the parameters are 
the following: d ~ 10 -f- 100 kHz, S ~ 100 ~ 500 kHz and 
A ~ 1 -J- 10 GHz. We propose to adiabatically eliminate the 
upper manifold and to study the dynamics of the lowest one 
with an effective Hamiltonian H per t- 



submanifold: 



ilpert / ^ ^ 



JSial 2 , |5i, 2 | 2 



6-d 



\S 2 ~ 12 



\Si, 2 \ 2 



S-d 

2 



IS!-' 2 



5 -2d 
1,3)(1,1| + H.c 



[|1,1)(1,1| + |1,3)(1,3| 



1,2><1,2| 



Sl,2 Si A f 1 



1 



Adt 



2 \S-d 5 -2d 

^2.2 Slfi ( 1 1\ idt 

2 \S-d + Sj e 
|1,2)(1,1| + |1,2)(1,3|] + H.c. 



(A4) 



Using this Hamiltonian we study the transfer rate of pop- 
ulation from level |1, 1) and |1,3) and viceversa. The 
main contribution is the direct coupling 



J 



(i) 

13 



|Sl,2| 2 ^ 2 

8 



(A5) 



A second contribution, which in our system will prove to 
be not-negligible, comes from a sort of "adiabatic elimi- 
nation" of the level |1,2): 



7(2) _ _ 

J 13 ~~ 



{l,3\H pert \l,2) (1,2\H 2 



pert 



i,D 



(l,2\H„ 



pert 



1,2) 



(l.llflpertll.l) 



(A6) 



The main approximations in the presented "6-levels 
model" reside in the neglection of delocalized higher- 
energy free states and in the neglection of couplings be- 
tween neighbouring higher-bands localised wannier func- 
tions. These processes could both induce spurious pop- 
ulation transfers to next-neighbouring states. Regarding 
the first issue, this mainly means that experimentally 
there is a trade-off between a large detuning regime, al- 
lowing powerful lasers and strong effective couplings with 




Population Level ll,l> 
Population Level 1 1 ,3> 



20 
time [ms] 



40 



FIG. 12: Exact time evolution of the populations with F z = 
of the 6 x 3-levels model describing the dynamics of the 
hyperfine levels of the ground state under the action of three 
Raman couplings. The initial state is |1, 1). Inset contains the 
maximum population reached in each level. The parameters 
used are listed in Table HV] 



noisy spurious population transfers, and a small detuning 
one, with small clean couplings. Concerning the second 
point, this essentially implies a trade-off between a deep 
lattice configuration with localized wavefunctions and a 
shallow lattice regime, with spread ones. In the former 
case neighbouring wannier functions are not connected by 
the Hamiltonian, but the overlap between different bands 
Si t 2 is also negligible; in latter case the Si 2 becomes im- 
portant but couplings between neighbours become also 
significative. 

It is possible to engineer the previous setup in order 
to get an independent control on the hopping rates of 
the different spin species. To this aim, we must be able 
to independently couple desired pairs of states (physi- 
cal F = 1 and auxiliary F — 2 states) via independent 
Raman transitions. This can be realized with the help 
of "energy selection rules" , i.e. choosing pairs of states 
with distinct energy differences and with the help of Ra- 
man transitions far-detuned from all the energy differ- 
ences but that of the pair that they should couple. In 
our case we take advantage of the fact that the hyper- 
fine Lande factors of the manifolds F = 1 and F = 2 
are one the opposite of the other, and suggest to split 
the hyperfine manifolds with a magnetic field and to use 
the level \a = (F = 2, F z = m)) as ancilla state for the 
\a = (1,to)) level. The energy differences of these pairs 
can then be detuned of circa 100 MHz with moderate 
fields of circa 66 G; in this case the use of Raman cou- 
plings detuned from the pair transition of hundreds of 
kHz or even MHz would do the job (see Fig. El). 

We now report on some numerical simulations which 
consolidate the arguments given above and show that it 
is indeed possible to induce different hopping rates for the 
three spin species. We calculate the exact time evolution 
of a 6 x 3-levels model under the action of three differ- 
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FIG. 13: Time evolution of the |1, 1) levels with different mag- 
netic numbers. As explained in the text, with the parameters 
listed in Table |TV] it is possible to tune the different hopping 
rates to very different values. 



ent Raman couplings. The parameters characterizing the 
levels and the couplings are listed in Table [TV| We assume 
the possibility of engineering Raman couplings carrying 
no angular momentum (w transitions) and therefore ne- 
glect the possibility of transferring population between 
states with different F z . This factorizes our 18-levels 
problem into three 6-levels problems, which are numer- 
ically less demanding. In case this were not possible, 
spin mixing would still be almost prevented by energy- 
conservation constraints. 



We show in Fig. 12 the exact time evolution of the six 
levels with F z = 0. At the beginning all the population is 



Level & F z 


Energy 


Parameters 




1,1), m 


fiFBm 


Ahf 


6.8 GHz 


|1,2), m 


^.pBra + 5 


fj, F B 


50 MHz 


|1, 3), m 


fipBm 


8 


60 kHz 


|2,l),m 


Ahf — fiFBm 


Si, 2 


0.2 


|2,2),m 


Ahf — fiFBm + S 


Si,i 


0.6 


|2,3), m 


Ahf — fiFBm 


S2,2 


0.6 



# Raman 


D. 


LU 




ang.mom. 


1 
2 
3 


31 kHz 
45 kHz 
39 kHz 


A H F-2fi F B + S- 
A H f + 5 - 
A H F + 2fi F B + S- 


- 300 kHz 

- 300 kHz 

- 300 kHz 








F z 



J 



(i) 



J 



(2) 
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-135.0 Hz 
-270.0 Hz 
-202.5 Hz 



-18.3 Hz 
-73.0 Hz 
-41.4 Hz 



Extimated T Numerical T 



20.4 ms 
9.1 ms 
12.8 ms 



22.1 ms 
10.1 ms 
14.0 ms 



TABLE IV: Parameters used in the simulation of the 3x6- 
levels model. We compare the periods of the Rabi oscillations 
with the theoretical values calculated taking into account only 
the Raman coupling which is supposed to drive the transi- 
tion. Discrepancies are contributions of the off-resonant Ra- 
man couplings and higher order corrections. 



in the level |1, 1) and very clear Rabi oscillations between 
the |1, 1) and |1, 3) levels can be seen. It is very important 
to notice that only a negligible fraction of the population 
is lost into the other four states (see the inset in Fig. 12 1. 
In Fig. [l3]we compare the time evolution of the levels 
\F Z = —1, 0, +1, k = 1): it is very interesting to observe 
that it is possible to induce different hopping rates for 
the three spin species. Indeed, as shown in Table IV the 



simple application of Eq. [A5| and [A6| corroborates the in- 
tuition that each hopping rate is ruled by only one Raman 
coupling, whereas the action of the others, far detuned, 
introduces only small corrections. It is also possible to 
check that including the other lasers the agreement with 
the experimental data improves. 

Concluding, this analysis shows that optical lattices 
loaded with alkaline atoms display a hyerarchy of energy 
scales which could be exploited to engineer hopping op- 
erators breaking the SU(2) symmetry. In particular, we 
propose to employ the upper hyperfine manifold to pro- 
vide auxialiary states to be adiabatically eliminated and 
to take advantage of superlattice configurations to trap 
them in the middle of each link. 



Appendix B: Magnetic Flux Quantization Condition 



In Sec. VI we deal with a (discretely) translational in- 



variant two dimensional lattice pierced by an external 
(homogeneous) magnetic field. In this Appendix we pro- 
vide more details on the study of such a system via a finite 
lattice with periodic boundary conditions (PBC). In par- 
ticular, we show that the need for mutual commuting 
Hamiltonian and discrete-translation operator imposes 
non-trivial conditions on the dimension of the sample. 



1. Bulk 

We start discussing the Hamiltonian and the discrete 
translation operator in the bulk. We consider the Landau 
gauge: A = B(0,x) and introduce the number of fluxes 
crossing each plaquette a = Ba x a y /&o where $0 is the 
flux quantum and a x and a y are the dimensional lattice 
constants. From now on x and y will just be adimcnsional 
integer numbers labelling the sites of the lattice. 

The standard generalization of the Bose-Hubbard 
Hamiltonian in presence of an external magnetic field is 
the Harper Hamiltonian: 



H 



-J 



-1"k%olx i\ 






h x,y+l u, x,y 



F d 



iJ.c.(Bl) 



where d x>y and d x y are boson annihilation and creation 



operators satisfying [d x . y ,d x , , 



Oxx'V' 



vv' ■ 



The action of the standard discrete translation opera- 
tor T mn — T{ma x + na y ) (m, n € N) on the field opera- 
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tors is the following: 



1,0 u x ,y J L0 " 


- d (t) 


0,1 u x,y J 0,l " 


- d (t) 


T 


mm Tin 
" J l,0 J 0,l 



J 0,l J l,i 



Since these T TO: „ operators do not commute with the 



Hamiltonian in Eq. Bl we need a "magnetic" translation 
operator M m>n commuting with the Hamiltonian, which 
are the discrete version of the continuum case discussed 
in Ref. WE: 



M 0) i 4" Ml, 

lvl m.n^'x,y lvl m,n 



+(-)2may >(t) 
e tt x+l,: 



d (t) 
u x,y+l 

e -(+^ amn M m M^d^yMll m\™ 
e +(-^ amn M n^ M™ d^ v M\™ Mly 



Last equation indicates clearly the peculiarity of the mag- 
netic translations which leads to the Aharonov-Bohm ef- 
fect: the result of a translation from one point to another 
one strongly depends on the followed path and eventu- 
ally, moving along a closed loop, gives to the state a phase 
proportional to the encircled magnetic flux. 

We can verify the commutativity of the magnetic 
translation operator with the Hamiltonian just by check- 
ing Miq because translations along other directions 
commute straightforwardly: 

M 1)0 e-™ a *dl >y+1 d x , y M 1)0 = 

— „-2iriax -2iria(y+l) +2-rriay jf J 

— c e e u x+ 1,8/4-1"'*+ 1,1/ 

_ -2TTia(x+l) A ■ 

— c "a;+l : j/+l u ^+l,i/ 

^i,o d x+ly d x , y M 10 = d x+ltV d XtV 
^ M lfi H M 1]0 = # 

where we exploit the sum over dummy x in H and change 
variables to x' = SE + 1, always possible in the bulk. 

2. Boundaries 

We now discuss the possibility of studying the pre- 
vious infinite system with a finite system of dimension 
L x x L y with PBC. L x . y are here adimensional numbers 
which can be used to define the total number of fluxes 
crossing the finite system: N$ = L x L y a. In order to be 



able to identify the bosonic operators residing on sites 
whose distance is mL x a x + nL y a y , with m,n £ N, we 
must require the total number of fluxes N$ to be an in- 
teger number. This can be proven simply translating 
the field operator around one plaquette L x x L y . As be- 
fore, we also require the Hamiltonian and the "magnetic" 
translation operators to commute; in particular, we dis- 
cuss in detail the interesting case of translation along x: 
Mi )0 H Ml Q = H. 

We separately analize this equation on each link of the 
finite lattice. In particular, when considering the links 
oriented along the y direction, it reduces to the following 
equality: 



— 2-Kiax -2iria [y+l-y] jf j _ 

e e ' a xTi.^+i X + 1 ^V ~ 

— p -2Triax+l jf J 

~ e a ^+i i y+i a x+l,y 



(B2) 



where x + 1 denotes the modulus count (x - 
the same holds for y. 

We distinguish four cases: 



1 mod L x ); 



1. x G [0, L x - 2] A ye [0, L y - 2]: Eq. JB2] is au- 
tomatically satistified, as it happens in the bulk; 



2. x = L x -1 A y G [0,L„-2]: Eq. B2 is fulfilled only 



if e -27ria(L a! -l) e -27ria _ ^ wh j ch imp l ies aLx g 



3. » G [0, i x - 2] A y = L y - 1: Eq. B2 is fulfilled only 



if e +2*-*a(£„-i) = e ~ 27rM , which implies aL„ G N; 



4. a; 



1 A 



.'/ 



only if e-^MLv-i 
a(L y - L 8 ) G N. 



2/ 

) „+27rza 



1: Eq. B2 is fulfilled 



( l m i) — i ; which implies 



The double constraint N<j>/L x , N$/L y G N and the 
desired magnetic filling one N$ — N strongly reduces 
the number and variety of finite size systems numerically 
treatable with moderate effort. The Hilbcrt space for 
the examined 4x4 lattice with 4 particles consists of 
3.620 states, but already 5 three-hardcore bosons on a 
5x5 grid need 110.630 to be described. Such strict con- 
straints could be circumvented if one introduces proper 
singularities of the magnetic field to fulfill the correct 
translational and periodic conditions; however, any spu- 
rious correction introduced by hand would strongly affect 
the numerics on the small scales treatable. Therefore we 
decided to stay stuck to the strictest version given above. 
An extensive numerical study of this problem, though in- 
teresting, goes well beyond the purposes of the present 
paper and is left for future investigations. 
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